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ABSTRACT 

We discuss the effects of inhomogeneous sky-coverage on CMB lens reconstruction, fo- 
cusing on application to the recently launched Planck satellite. We discuss the "mean- 
field" which is induced by noise inhomogeneities, as well as three approaches to lens 
reconstruction in this context: an optimal maximum-likelihood approach which is com- 
putationally expensive to evaluate, and two suboptimal approaches which are less in- 
tensive. The first of these is only sub-optimal at the five per-cent level for Planck, and 
the second prevents biasing due to uncertainties in the noise model. 

Key words: cosmic microwave background - methods: numerical - cosmology: ob- 
servations — gravitational lensing 



1 INTRODUCTION 

The current generation of Cosmic Microwave Background 
(CMB) data has proven to be remarkably well approximated 
as a statistically isotropic Gaussian random field (Komatsu 
|et al.|[2009[ ). Upcoming experiments are expected to push 
decisively past this approximation, to reveal a CMB which 
has been subtly distorted by gravitational lensing due to the 
large-scale-structure (LSS) which intercedes between our- 
selves and the surface of last scattering ( Lewis & Challinor 
2006). Mathematically, the effect of a fixed LSS realization 



is to make the CMB statistically anisotropic, introducing 
off-diagonal elements into its covariance. This enables one 



to construct estimators for the lensing potential ( Hirata fc 
|Seljak||2003| |Okamoto fc HupOOl} . The power spectrum of 
the measured lensing potential may then be used to obtain 
improved parameter constraints, particularly for parameters 
which affect the late-time evolution of the Universe. The re- 
cently launched Planck satellite, for example, is expected to 
measure the CMB lensing signal internally with cosmolog- 
ically useful precision, enabling it to constrain the sum of 
neutrino masses to O.leV (iLesgourgues et al. 20061. Lens- 



ing is also important as a potential contaminant for non- 
Gaussianity studies. The cross-correlation between the lens- 
ing potential and ISW/Rees-Sciama induced temperature 
fluctuations results in a bispectrum which has large overlap 
with the "local" type of non-Gaussianity. Lensing therefore 
results in a bias to primordial non-Gaussianity estimation, 
which will be significant for Planck ( Serra fc Cooray|2008 1. 
Correction for this bias will be aided by an accurate lensing 
reconstruction. 



CMB lensing reconstruction works on the assumption 
that the underlying CMB is statistically isotropic, and that 
any statistical anisotropy is due to gravitational lensing. 
As such, it is potentially contaminated by any systematic 
which introduces anisotropy onto the observed sky: beam 
asymmetries, astrophysical foregrounds, and inhomogeneous 
sky-coverage are all expected to complicate the lens recon- 
struction. In the absence of computationally expensive de- 
convolution mapmaking ( Armitagc & Wandclt 2004), beam 
asymmetries represent an unavoidable source of systematic 
error, which will need to be quantified for any ultimate lens- 
ing analysis with Planck. Foregrounds, on the other hand, 
may be cleaned to a high degree of accuracy by exploiting 
Planck's wide frequency coverage. The magnitude of resid- 
ual foregrounds at the small scales of interest to lensing re- 
construction will not be adequately understood until Planck 
has started to collect data, however. In this work, we will fo- 
cus on the effects of inhomogeneous sky-coverage. The scan 
strategy of Planck will result in noise levels which depend 
strongly on ecliptic latitude, and so the effects of noise in- 
homogeneities are a large concern. We will study both the 
optimal treatment of noise inhomogeneities, as well as two 
suboptimal approaches: one which is computationally sim- 
pler than the optimal reconstruction, and one which is in- 
sensitive to the instrumental noise model. 
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1.1 Lens Reconstruction 

We begin by reviewing the methodology of lens reconstruc- 
tion. Consider a data model given by 



6(0) = 6(0 + V$(«)) + n(Q) 



(1) 



where is the observed CMB, O picks out a loca- 
tion on the unit sphere, 3> is the lensing potential (for 
more details see Lewis & Challinor 2006), B is the pri- 
mary, unlensed CMB with covariance 7P e , and n(Q) is 
the instrumental noise realization with covariance ma- 
trix C nn . Throughout this work, we will use a fixed 
flat ACDM cosmology for C ee , with standard parame- 
ters {n b ,Cl c ,h,n 3 ,T,A s } = {0.05,0.23,0.7,0.96,0.08, 2.4 x 
10" 9 }, wh ich is consistent w ith the WMAP5 best-fit power 



spectrum ( Nolta et al.|2009 l. The maximum-likelihood esti- 
mator for the CMB lensing potential in the limit of small "I> 
is given by 



4>z 



(2) 



where A is a normalization matrix, the average is taken over 
realizations of the CMB and noise, and the un-normalized 
estimator $ is given in harmonic space by 

I I' 
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Here 6 = (C ®)' 1 ® = (C ee + C nn )- 1 Q is the inverse- 
variance filtered sky-map. The indices L and M give the 
mode of the lensing potential which is being reconstructed. 
The geometric term f lL ii is given (with the notation E; = 
I 2 + I) by 
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The estimator normalization matrix A is also eq ual to the 
estimator covariance JV**. For more details, see Hirata & 



Scljak 2003s |Smith et al]2007| 

This likelihood motivated estimator is closely related 
to the minimum-variance quadratic estimator of Okamoto 



and Hu (Okamoto & Hu 20031, which may be derived un- 



der the assumption of uniform sky-coverage (homogeneous 
noise and no masking), in which case C nn is diagonal. The 
likelihood approach motivates two modifications which im- 
prove the performance of the estimator for non-uniform sky- 
coverage: 

• Mean field subtraction: The effect of non-uniform sky 
coverage is to introduce off-diagonal elements into the har- 
monic space noise covariance matrix which are interpreted 
by the estimator as lensing effects and give the estimator a 
non-zero expectation even in the absence of lensing. 

• Anisotropic inverse-variance filtering: This is an intu- 
itive generalization from the rotationally invariant filters of 
the Okamoto and Hu estimator, although difficult to derive 
in that context. 

Similar ingredients are seen e.g. in bispectrum estimation 
(Creminelli et al.||2006"|. Application of (C ee ) _1 is gen- 



et al.|[2007 l, however with full-sky coverage we find that it 
is sufficiently well conditioned that it may be applied using 
conjugate-gradient descent with a diagonal preconditioner 
in less than one hour on a 2GHz processor to £ max = 2500, 
with an average fractional error of O(10~ 6 ) for each mode 
of the inverse- variance filtered field. 

As a baseline, we will also consider approximating the 
inverse variance filter as rotationally invariant, taking only 
its diagonal elements, averaged over the azimuthal index m. 
In this case, we take C® § = C ee + diag(C nn ), where the 
diagonal operation is given by 



diag(C7) i 



^ll'^mm' 2^ _j_ ^ ^ ^lm" ,lr, 



(5) 



We will refer to this as the "uniform" estimator, and symbol- 
ically denote it in lowercase as cf>. The corresponding normal- 
ization and covariance matrices will accordingly be denoted 
as A and respectively. Expressions which involve cj) and 
A will implicitly be taken to use the symmetrized inverse 
variance filter as well. The normalization to lens fluctua- 
tions is diagonal, independent of M, and may be calculated 
analytically ( jOkamoto fc Hu|2003[ ): 



At = (2L+ I) 
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(6) 



erally a challenging problem at Planck resolution (Smith 



With this normalization and accurate mean-field subtrac- 
tion, the uniform estimator produces an unbiased recon- 
struction of the CMB, however the normalization and the 
estimator variance are no longer explicitly equal (although 
we will see that in practice they are still very close). 

1.2 Noise Model 

To illustrate our discussion of inhomogeneous noise effects 
on lens reconstruction we will work with a semi-realistic 
model for Planck using simulated data from the detectors at 
143GHz, with an isotropic Gaussian beam of ofwhm = 7', 
at HEALPix N side = 2048. 

We will assume that the map noise is Gaussian and ef- 
fectively uncorrelated between pixels (prior to beam decon- 
volution). In reality, non- white noise below the instrumental 
1// knee frequency leads to inter-pixel noise correlations in 
the Planck data. To accurately study these effects in the con- 
text of lens reconstruction requires many simulations gen- 
erated by performing the mapmaking procedure on realistic 
time-ordered data, a computationally expensive task. We 
leave this study to a future work, concentrating for now on 
pixel-uncorrelated noise. In this case, the noise is completely 
characterized by a variance map. We obtain this map from 
the output of the Springtide destriper mapmaker applied to 
a single realistic Planck simulation ( jAshdown et al.|007a|b[ ). 
The noise levels which result have a white power spectrum 
which is twenty per-cent greater than that for the "Blue- 
book" value of 43/i7f • arcmin ( [Efstathiou et aL]|2006[ ), re- 
sulting in a cosmic variance limit of approximately £ = 1500. 
This enhancement of the white-noise noise level is due to the 
inhomogeneity of the sky coverage. 

The noise variance map is displayed in the upper panel 
of Fig.[T] The features in the noise variance are due primarily 
to the Planck scan strategy. From the L2 point of the Earth- 
Sun system, Planck spins at approximately one rotation per 
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minute about the anti-Sun direction. The angle between the 
optical axis and the spin axis is 85°, and so the detectors 
trace rings which are nearly great circles on the sky. This 
results in relatively low noise levels near the ecliptic poles, 
where there are many observations as well as a large degree 
of cross-linking. To gain sky coverage at the poles them- 
selves, the Planck model which we have used employs the 
so-called "cycloidal" scan strategy, in which the spin axis 
inscribes a circular path around the anti-Sun direction with 
a period of six months, keeping the angle between the spin 
axis and the anti-Sun direction at 7.5°. This results in the 
cusp features near the ecliptic poles in the upper panel of 
Fig. [I] The thin circle of low-noise levels which connects the 
ecliptic poles is because the simulation which was used as 
input for the mapmaking procedure included just over one 
year of data, and so this section of the sky has on average 
fifty per-cent more hits than other regions. 



2 RESULTS 
2.1 Mean field 

We begin by considering the mean field term. For the uni- 
form estimator, the corresponding mean field may be calcu- 
lated analytically: 
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This expression may be reduced using the Gaunt integral 
and the orthogonality properties of the Wigner-3j symbols 
to give a simplified expression for the mean field: 
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,(8) 



where Nlm is the harmonic transform of the noise-variance 
map, n p i x is the number of map pixels and Bi is the instru- 
mental beam transfer function. We can see that for the uni- 
form estimator, the mean field is simply the noise variance 
map convolved with a rotationally invariant response filter. 
In Fig. [T]we plot maps of the mean field for our Planck noise 
model. In Fig. [2] we plot the power spectrum of the mean 
field, as well the expected lensing power spectrum and es- 
timator variance. At low multipoles the magnitude of the 
mean field is considerably larger than the estimator vari- 
ance for uniform noise. 

For the anisotropic estimator, the mean field term is 
even larger. In addition to the noise-only mean field, the 
anisotropic filtering generates a mean field from the CMB 
anisotropics themselves. We will now consider the variance 
of the lensing reconstruction after subtraction of this large 
mean. 



2.2 Estimator variance 

For homogeneous white noise with a power spectrum given 
by Eq. jH]), the uniform estimator variance is equal to its 




Figure 1. Noise variance map (top) and lens reconstruction mean 
field (lower two panels), for our Planck noise model, in galactic 
coordinates. The thin stripe of low noise levels through the ecliptic 
plane is because we have used a cosmology in which data formats 
dominate over aesthetics at z = 0. 




Figure 2. The power spectrum of the mean field due to inho- 
mogeneous noise for the uniform estimator (solid), the fiducial 
estimator variance (dashed), and the expected cosmological 

power spectrum of the lensing potential (dotted). 
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normalization, given in Eq. (fsj) . We will therefore refer to 
this as the fiducial estimator noise, Nt£ = A. 

In Fig. [3] we compare the variance of the mean- 
subtracted uniform estimator with this fiducial value. They 
agree well within our Monte-Carlo error bars for £ < 300. 
The explanation for this agreement comes from considering 
the anisotropic, Gaussian noise distribution instead as an 
isotropic, non-Gaussian field, e.g. by randomizing the orien- 
tation of the scan strategy. For the uniform estimator, the 
isotropic component of the variance (equal to N££ ) is unaf- 
fected after mean-field subtraction, however the anisotropy 
of the noise distribution before averaging over orientation re- 
sults in non-Gaussian connected terms. The "primary" con- 
figuration of this trispectrum is removed by subtraction of 
the noise mean field, and at low-£ the remaining configura- 
tions are suppressed by the estimator, as noted in |Okamoto| 
Ifc Hull2002l and IHanson et al. 120091 

For the anisotropic estimator, it is generally most effi- 
cient to compute the estimator covariance by Monte-Carlo, 
using 



Cov[3>- ($>] = (iV**)~ 



(9) 



keeping in mind that in our notation, the $ estimates are un- 
normalized. For a small number of parameters, estimation 
of this covariance matrix and subsequent inversion is typ- 
ically stable. In the bispectrum context, for example, one 
usually seeks to determine the projection of a prescribed 
bispectrum shape in the data, a single parameter. In the 
case of lens reconstruction, however, the Fisher matrix will 
contain thousands of useful modes and is effectively impossi- 
ble to obtain with sufficient accuracy from Monte-Carlo. In 
practice, this may not be an issue as many uses of the recon- 
structed lensing potential require (A r **) _1< i > rather than $ 



itself (e.g. Smith et al.|2007[ ). In this work, we simply place 
a lower limit on the variance, using the result that 



diag(iV**) > [drngdN* 



(10) 



which holds for any covariance matrix. This estimate of JV** 
is also compared to the fiducial variance in Fig. [3] In prac- 
tice the correlations between the reconstructed modes, which 
have been completely neglected here, reduce the amount of 
information in the estimator. The optimal anisotropic filter- 
ing results in a lensing estimator with approximately five 
per-cent less variance than the uniform estimator, a poten- 
tially useful improvement, although it must be kept in mind 
that this is an optimistic result. 

At low multipoles we have seen that the magnitude of 
the mean field is considerably larger than the estimator vari- 
ance. If the noise is well-understood, however, then the in- 
crease in the variance of the mean-subtracted estimates is 
negligible, for both the optimal anisotropic and uniform fil- 
tering approaches. As a possible systematic check, in the 
next section we consider a lensing estimator which is insen- 
sitive to noise inhomogeneities. 
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Figure 3. Variance of lensing reconstruction with three differ- 
ent estimators, compared to the fiducial noise level: the "uni- 
form" estimator (circles), the lower-limit of Eq. | |10| l for the opti- 
mal "anisotropic" inverse-variance filter estimator (squares) and 
a pessimistic result for cross-maps estimator (diamonds). Each 
point is the average for 100 Montc-Carlo simulations. 



with independent noise the "optimal" treatment is to con- 
dense them to an inverse variance weighted average, which is 



then analysed as a single map (see e.g. Hamimeche & Lewis 
2008). To avoid biasing due to uncertainties in the noise 
model, it can be useful instead to work strictly with pairs 
of maps, such that auto-correlations of the noise are never 
produced. In the context of lens reconstruction, this allows 
us to avoid the mean field due to noise anisotropies. Given 



two noisy maps ai r , 



we consider the cross-estimator 
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where the Di are isotropic filters and the normalization to 
lens fluctuations is given by 



A a L xb = (2L + 1) 
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This type of estimator has been used by Hirata et al. ( 2008 \ 
for the purpose of cross-correlating <f> with large-scale struc- 
ture from galaxy surveys. Here we are more interested in the 
lensing potential power spectrum, which can be estimated 
as 



c = 



2L + 
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2.3 Cross-maps 

In this section we consider the behaviour of the lensing esti- 
mator on sets of maps with independent noise realizations. 
From the likelihood approach, given a set of CMB maps 



where 5 is a collection of map quadruplets and \ S\ is its size. 
For simplicity we will assume that all of the maps have the 
same noise properties (although different noise realizations) 
and so the D; filters should all be equal. In this case the re- 
construction noise bias to C"*'*' for homogeneous noise would 
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be given by 



3 CONCLUSIONS 
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where A L is Eq. ( 12 1 evaluated for the common Di fil- 



ter. Noise anisotropies will manifest themselves in additional 
contributions to N x , however we have already seen in Sec- 
tion [22] that these terms are small at low-^. The choice of 
S determines how sensitive this estimator is to instrumental 
noise and inhomogeneities thereof. 

If we have n independent maps and S is taken to con- 
tain the iC n quadruplets for which (a ^ b ^ c 7^ d), then 
it can be shown that N xx is minimized for Di = Cf e 
and is equal to the reconstruction noise level for a cosmic- 
variance limited experiment. This is directly analogous to 
the cross-correlation approach which is often used for tra- 
ditional power spectrum estimation, as it removes any in- 
strumental noise bias from the Cf estimates. The low re- 
construction variance is somewhat misleading however, as 
the variance of the power estimates contains contributions 
from the instrumental noise which must be accounted for 
in a parameter analysis. This approach discards a fraction 
n!/[(n — 4)!n 4 ] of the possible map combinations, and so its 
effective sensitivity must be less than the optimal approach 
of performing reconstruction on the minimum-variance sum 
of the maps, particularly for small n. 

If we take S to be the set of all quadruplets with 
(a^=b)(c^d) then we retain more of the possible map com- 
binations. Cf/ is no longer completely free of noise depen- 
dence, however the noise only enters through the map 
power spectra, which are experimental "observables" and do 
not rely on any modelling of the noise. In Fig. [3] we plot the 
variance of the cross-estimator measured from simulations, 
relative to the minimum-variance result for the pessimistic 
case of only two maps. We generate realizations of ai m and 
bi m with twice the noise variance levels of the previous sim- 
ulations, such that the noise level of the minimum-variance 
average is unchanged. We inverse-variance filter them with 
Di = Cf 8 + 2C nn , which minimizes the term with the 
largest contribution to iVjJ, . If the D filter were permit- 
ted to couple 1,1' then an estimator with smaller variance 
could be derived, however it does not have a known fast po- 
sition space form to make its calculation feasible at Planck 
resolution, and so we do not consider it here. In any case, 
our purpose in this section is to demonstrate a consistency 
test rather than a minimum-variance reconstruction of 4>. 
The sub-optimality of this approach is evident, with the 
variance of these estimates being approximately twenty per- 
cent larger than the fiducial value (although this discrep- 
ancy will be less for a larger number of maps) . The need to 
perform mean field subtraction has been obviated, however. 
The effect of the inhomogeneous noise is negligible at the 
level of the estimator variance in this approach. Only the 
noise and CMB power spectrum are required. Noise inho- 
mogeneities do make contributions to the higher moments 
of the reconstruction statistics, however, as can be seen from 
the increased size of the Monte-Carlo error bars. 



We have studied the effects of inhomogeneous instrumental 
noise on CMB lens reconstruction. The main effect is to 
introduce a mean field even in the absence of lensing, with 
large power on scales I < 100. 

We have studied the optimal estimator in the case of 
such inhomogeneities, and found that for Planck it performs 
approximately five per-cent better than a suboptimal ap- 
proach which is computationally less expensive and easier 
to study analytically. With accurate mean field subtraction, 
the suboptimal "uniform" estimator itself performs as well 
as would be expected for homogeneous noise with the same 
power spectrum. 

Both the optimal and uniform reconstructions require 
an accurate modeling of the noise inhomogeneities to be ef- 
fective. This requirement may be bypassed using a lens re- 
construction based on pairs of maps with uncorrelated noise, 
which we suggest will provide a useful consistency test. 
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